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Abstract 

The VETA-I mirror was calibrated with the use of a 
collimated soft X-ray source produced by electron 
bombardment of various anode materials. The FYVHM, 
effective area and encircled energy were measured with 
the use of proportional counters that were scanned with 
a set of circular apertures. The pulsers from the 
proportional counters were sent through a multichannel 
analyzer that produced a pulse height spectrum. In 
order to characterize the properties of the mirror at 
different discrete photon energies one desires to extract 
from the pulse height distribution only those photons 
that originated from the characteristic line emission of 
the X-ray target source. 

We have developed a code that fits a modeled 
spectrum to the observed X-ray data, extracts the counts 
that originated from the line emission, and estimates the 
error in these counts. The function that is fitted to the X- 
ray spectra includes a Prescott function for the resolution 
of the detector a second Prescott function for a pileup 
peak and a X-ray continuum function. The continuum 
component is determined by calculating the absorption 
of the target Bremsstrahlung through various filters, 
correcting for the reflectivity of the mirror and 
convolving with the detector response. 


Introduction 

The VETA-I (Verification Engineering Test Article-I) 
consists of a pair of Wolter Type-I mirrors which are 
intended to be used as the outermost set of the AXAF 
(Advanced X-ray Astrophysical Facility) telescope. A 
detailed description of the performance and scientific 
advances anticipated from AXAF observations are listed 
in reference 1 . The improved sensitivity of the AXAF 
mirrors compared with those of the Einstein and ROSAT 
X-ray telescopes results mainly from the larger 
attainable effective area of AXAF over a broader energy 
bandwidth and the higher angular resolution resulting 


from improvements in the mirror figure and surface 
smoothness. A natural consequence of these improve- 
cents is the need to use detectors with higher spatial and 
energy resolution and the need to employ more 
sophisticated data reduction techniques. The physical 
properties of the VETA-I optic measured are the FWHM, 
the effective area and encircled energy. Definitions of 
these quantities are given in Kellogg et al 2 . For the initial 
analysis of the VETA-I proportional counter spectra we 
used a quick-look analysis technique which basically 
involved the summation of counts in a selected region of 
interest. The quick-look data provided a fair description 
of most of the mirrors properties with a moderate 
amount of analysis. In this paper we will describe a more 
elaborate method for analyzing the VETA-I data which 
leads to a determination of the mirror effective area to 
better than 5% and the encircled energy better than 2%. In 
our analysis we will apply a set of corrections for the 
spectral contamination of the measured proportional 
counter data due to the existence of a bremsstrahlung 
component in the X-ray source, pulse pileup effects 
which distort the spectrum, background and deadtime. 

In section 1 we provide a brief outline of the 
experimental setup. A detailed description of the model 
employed for the data reduction is given in section 2. 
Finally section 3 is devoted to the presentation of the 
results of the VETA-I data reduction and a discussion of 
the effect of these corrections on the measured mirror 
properties. 


1. 0 The VETA-I calibration set up 

The VETA-I calibration set-up is shown in figure 1. An 
electron impact X-ray source with interchangeable 
targets and filters is located 528 m away from the optic in 
order to approximate an incident plane wave. The 
parabola PI and hyperbola HI mirrors are made out of 
zerodur and were supported by flexures. Their relative 
alignment and centering were adjusted with a set of 
linear actuators. The X-ray flux incident on the VETA-I 
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was monitored by sealed and flow proportional 
counters manufactured by LND Inc. while an identical 
pair of sealed and flow counters were located at the focal 
plane of the VETA approximately 10 meters behind the 
optic. The flow counters have a graphite coated 1.7 Jim 
thick polypropylene window. For energies below 1.49 
keV the flow counters were filled with methane at a 
pressure of 250 Torr while for energies above 1.49 keV a 
mixture of 90% Ar and 10% methane (P10) at a pressure 
of 400 Torr was used. The electronics associated in 
acquiring the data are shown in figure 2. A set of 
periodic pulsers (not shown in the figure) were injected 
into 2 multichannel analyzer (MCA) channels in order to 
monitor the deadtime of the pulse height system. The 
focused X-ray image was scanned by translating a plate 
containing a set of pinholes, slits and annuli across the 
image with a position accuracy of approximately 1 4m. 

2,0 Model of the VETA-I monitor and focal plane 
proportional counter spectra 

The description of our model will trace the 
sequence of events traced by X-ray photons as they are 
produced through the deceleration of electrons in vari- 
ous targets, transmitted through filters to minimize the 
bremsstrahlung component, reflected by the PI -HI mir- 
rors and finally absorbed inside a proportional counter. 
The resulting pulses are fed into a multichannel ana- 
lyzer which converts their relative pulse height to a 
channel address. 


2*1 X-ray source spectrum 

For the analysis of the VETA-I data we have 
adopted Kramers relationship for the bremsstrahlung 
emission of the soft X-ray source: 

1(E) = kZ(E max -E) (1) 

where E max is the high voltage of the source anode, Z is 
the atomic number of the target and k is a fitted param- 
eter with units of inverse energy. The spectra that we 
obtained from the beam monitor detectors allow a 
direct fit to the continuum component at energies above 
the characteristic X-ray lines. 

22 Transmission of source filters and counter 
windows 

In section 2.6 we will show that one of the main con- 
tributors to the error in the derived focal plane (XDA) 
and beam monitor detector (BND) events is the contin- 
uum component under the line peak. In order to mini- 
mize the bremsstrahlung intensity under the peak a 
source filter is used with an absorption edge slightly 
higher in energy than the characteristic line. The prop- 
erties of the source filters used in the VETA-I test are 
shown in table (1). The errors quoted by the manufac- 
tures are near 30% and are due to errors in the measur- 
ing technique. The effective source filter thickness is the 
thickness determined by fitting a ensemble of BND 
spectra and performing a grid search of filter thickness 

vs. the ensemble average % as shown in figure 3. The 
68.3% confidence level is just x 2 min (unreduced) +1. 
Similarly an ensemble of BND spectra are fitted in order 
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Figure 2. Proportional counter electronics 
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to determine the strength of line and continuum emis 
sion of the X-ray source. 



The ratio of the effective thickness to the nominal 
value for the thicknesses ranged between 0.9 for Zr and 
approximately 2.6 for Mo. An independent measure- 
ment of the Mo thickness performed in a post VETA test 
provided similar results for the effective filter thickness 
and strongly suggests an error in the labeled Mo nomi- 
nal filter thickness. 


Table 1: 


Source 

filters 

Nominal 
Thickness pd n 

Effective filter thickness 
p d «ff 

Parylene 

2.20 x lO^gr cm’^ 

(1.7 ± 0.1) x ltT’grcm'* 

Ai 

2.70 x 10' 3 gr cm' 2 

(2.8 ± 0.1 ) x 10‘ 3 gr cm' 2 

Mo 

2.04 x lO^gr cm' 2 

(5.36 ± 0.05) x 10* 3 gr cm * 

Zr 

1.29 x 10 J gr cm ' 2 

(1.17 ± 0.03) x 10' 3 gr cm 2 

Cu 

4.48 x lO^gr cm* 2 

(5.82 + 0.5) x l(T*gr cm * 




2 

Figure 3. X grid search for source filter thickness , line 
strength and continuum strength 


23 Resolution, efficiency and uniformity of 
proportional counters used in the VETA-I calibration 


The resolution function R(E) of a proportional 
counter provides the response of a counter for an mo 
dent monoeneigetic beam of X-rays of energy E. Vanou. 
resolution functions have been determined theoretic* I iv 
by studying the statistics of secondary electrons pro 
duced by the primary photoelectrons and the multipi* 
electron avalanches initiated by each secondary electron 
For energies above 2 keV the response can be approo 
mated with a Gaussian with a mean proportional to th*- 

energy and a FWHM proportional to E 1/2 . This appron 
mation breaks down however at low energies where on*- 
observes an asymmetry in the response. This is mainlv 
due to the small number n of avalanches produced per 
secondary electron (n = the energy of the electron / ion- 
ization energy) and consequently for low energies one 
does not expect, according to the central limit theorem, 
to obtain a Gaussian distribution. 

For the present analysis the resolution of the propor- 
tional counter is approximated by a function first 

derived by Prescott 3 which takes into account the statis- 
tics of multiple avalanches. The Prescott function has the 
form: 


F(x) = 
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where the mean x of the distribution is x = m/Q and 

the standard deviation is G = Jim Q , Q is the Prescott 






width parameter. By minimizing the Prescott function 
F(x), -j^F(x) = 0, one can derive the following rela- 
tionship between the mean and the width parameter: 


m = .3Q 9 Q 2 

x 2 x 16 y 2 

max max X 


(3) 


In our algorithm the m and Q parameters of the line 
component of the spectrum are fitted parameters. For 
the determination of the measured continuum compo- 
nent one requires to know the width parameter vs. 
energy in order to convolve the detector response with 
the modeled incident continuum component over the 
measured energy band. 

By fitting a range of X-ray lines we have produced 
an empirical expression for Q(E) for E in the range of 0 
to 2.5 keV: 

Q ( E ) = a 0 E + a x E 2 + a 2 E? + a 3 £ 4 (4) 


with a 0 = 0.08514, = -0.13751, a 2 = 0.08026, a 3 =- 

0.01552 

The measured energy resolution and Prescott Q 
parameter for the proportional counters used in the 
VETA test are shown in figure 4. 


rates. Each pileup count corresponds to 2 events that 
were not separately detected in the multichannel ana- 
lyzer. At energies below the 1.49 keV A] line the poor 
detector energy resolution does not allow a direct deter- 
mination of the pileup component. The pileup peak rep- 
resents a small distortion to the spectrum and is cor- 
rected for by performing a simultaneous fit to the line 
and pileup peak. This fit provides a relation between the 
percent pileup counts and the total event rate. The mea- 
sured pileup versus total event rate is shown in figure 5. 
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Figure 5 Measured pileup rate vs. total MCA event rate 




Figure 4. Measured Prescott Q parameter and proportional 
counter energy resolution 


We have compared our results to a theoretical model 

presented in Johns and Yaffe 4 which predicts that the 
percent pileup follows the equation: 


N, 


N< 


^ = ( 1 . - e -2pT ) 


XDA 


(5) 


where Nxda are tota * number of focal plane events 
originating from the X-ray line, N pu are the photon 
events that are piled up, p is the average photon event 
rate (from all MCA channels) and % is the characteristic 
time of the MCA. The theoretical prediction is also 
shown in figure 5. The fitted value for t is 1.85e-7 sec. 
The measured pileup fraction f(p) is fit to a 3'rd degree 
polynomial: 

/(p) = b(\)p + b(2)p 2 + b(3)p 3 

where b(l) = 4.845e-7, b(2> = 5.876e-ll, b(3) = -2.441e-15 


2.4 Peak pileup correction 


2.5 Extraction of X-ray line counts 


Soft X-ray spectra obtained during the VETA -I cali- 
brations show evidence of peak pileup at high counting 


The function that describes the model fitted to the 
focal plane spectrum has the form: 



fxDA^ = fline( 0 +/<:*«/(') +/p«(0 ( 6 > 

where fxDA^ are the counts in the MCA pulse height 
channel i, f^e is the X-ray line component which 
includes all excited K and L lines and their escape peaks 
and is attenuated by the mirror reflectivity, f amt is the 
continuum part of the spectrum which also contains the 
fluorescence radiation that escapes the counter and the 
effect of the mirror reflectivity, f pu is the pileup compo- 
nent and f bgd is the normalized for live time back- 
ground. In the calculation of the continuum and line 
component we perform a convolution of the detector 
response R(Ej,Ej) with the incident X-ray spectrum. The 
incident spectrum is a sum of nl X-ray lines of strength 
a-^sfEj), where s(Ej) is the relative intensity of line Ej 

obtained from reference 3 , and a bremsstrahlung emis- 
sion BremfEfc) of strength a 2 - For the focal plane spec- 
trum we include the effect of reflectance off the mirror 
by multiplying with the model effective area A(Ej) as 
determined by ray tracing the PI -HI mirror surface and 
using Henke 90 optical constants for zerodur. The input 
model effective area of the VETA-I is shown in figure 6 
together with the energies at which the calibration was 
performed. 



Figure 6. Input model effective area of the VETA-I 


The line component of equation (6) is calculated by 
using the expression: 

nl 

ftuJf) = X (£,) 

j- 1 


x[ (1 -esc(E Jt E l ))R(E i ,E.) + 

+ esc (Ej, E,) /?(£,, E. - £,) J 

For the continuum component we used: 

ph 

X T^E k )T w (E k )Abs (E k )Brem(E k ) 

k= 1 

x[ (1 -esc(E k ,E,))R(E i ,E k ) + 

+ esc (E k , E ,) /?(£,, E k - E,) ] 


T f and T w are the transmission of the source filter and 
counter window respectively, Abs(Ej) is the efficiency of 
the counter, esc(Ej,Ej) is the probability of a fluorescence 
photon of energy Ej escaping from the counter times the 
fluorescence yield. 

The background component is normalized to take 
into account differences in the live time between the 
background run and the test run. f is the mea- 

sured background MCA spectrum. 


fbgd^ 


W-dt b gd) 

l bgd ( * - dt MCA l) 


'/ bgclW 


were t int is the integration time, t bgd is the integration 
time of the background run, dt bgd and dt^CA are the 
deadtime of the background run and test run respec- 
tively. 

For extracting the counts that originated from the 
X-ray line emission we begin by fitting the model func- 
tion of equation (6) to the spectrum using a Marquardt 
nonlinear least squares routine. We select a region of 
interest (ROI), determined by the XDA MCA channels 
(mo/mD and BND MCA channels (n^nD, around the line 
peak and sum the measured pulse height counts. We 
then subtract the fitted continuum component in the 
ROI, subtract the pileup in the ROI, add the fitted line 
counts outside the ROI, add two times the total pileup 
counts (every pileup count corresponds to two events 
detected as one count) and finally subtract the normal- 
ized background counts. The extracted line counts for 
the focal plane detector XDA p c b , corrected for pileup, 
continuum and background effects, are given by the 
expression: 

r - ■■■-■-- _ i 

171 j 171 q 1 

m V.*= X wo+ X wo+ 

i = m 0 i = 0 


( 7 ) 
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The extracted line counts for the beam monitor 


detector BND C d , corrected 

for continuum and 

back- 

ground are: 

n i 

n 0 - 1 


BND Cid = 

X $BNdW + 

X + 
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« =«o 

/ = o 



+ x fujA-i ui- 

i = + 1 i = n 0 

t inl X ( 1 -^) . 

(1- ~ dt BND ) bgt 

where Sxda^ and Sbnd( 0 are the counts in pulse 
height channel i from the total focal plane and beam 
monitor proportional counter spectra, dtxDA anc * dtgND 
are the deadtime in the MCA channels that contain the 
XDA and BND spectra respectively. 

The BND count rate is significantly lower than the 
XDA count rate and pileup effects are less than 0.1% 
and for this reason we have not included a pileup cor- 
rection for the beam monitor spectra. 


2.6 Error analysis implemented in spectral extraction 

The objective of the VETA-I calibration was to mea- 
sure the encircled energy with 2% precision and the 
effective area to 5%. In the present analysis we will 
present the calculated 1 a error components of the mea- 
sured effective area that are due to fitting errors, pileup 
effects, bremsstrahlung, continuum, background sub- 
traction and the deadtime correction. We now proceed 
in propagating errors through equations (7) and (8) that 
give the focal plane and beam monitor line events, in 
terms of the corresponding quantities in counts and 
percent deadtime. 


Subscripts d,p,c and b to a quantity refer to dead- 
time, pileup, continuum and background corrections to 
that quantity. XDA p c b are the counts in the XDA spec- 
trum originating from line emission and corrected for 
pileup, continuum and background. The error contribu- 
tions from each effect can be seen in the following equa- 
tions where we have neglected errors in the dead time 
for the background, which contributes less than 0.05%. 
The error in the focal plane counts is: 


= X X 
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For the error in the beam monitor counts we neglect 
the pileup terms due to the relatively low count rates 
measured in the BND detector. 
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a pu is the standard deviation in the pileup counts, a pum 
is the standard deviation of the pileup counts in the 
region of interest (m^m]). The errors due to the dead- 
time correction are: 



X $XDaW X ^ XDA O') 

' — /n^ i — wij 

t inl X ^ Pi 


°/pl + 


source filter thickness, line strength, continuum 
strength, the Q Prescott parameter of the detector 
response, the m Prescott parameter and the model 
scaling factor. 


3.0 Results and discussion 
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The pulse height channels over which the puls- 
er counts appear in the XDA and BND MCA spec- 
tra are (morris) and ( 114 , 1 x 5 ) respectively. Since the 
pulser used for the deadtime correction was peri- 
odic, the error in the pulser counts as measured in 
the MCA is : 


f(jt-4) 2 p (x )dx 

1 o 


0.289 counts 


while the error in the injected pulser rate is : f p * 
0.289/tpuiser/ where tpu^ is the time over which 
pulser rate is measured. 

The errors in determining the continuum com- 
ponent from the fit are: 
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The summations are performed over the np vari- 
ables (a(j) for the XDA spectrum and b(j) for the 
BND spectrum) of the continuum component; 


In figure 7 we present fits and their residuals to 
the beam monitor and focal plane spectra for Al, 
Zr,C and Mo source targets. The continuum com- 
ponent for the carbon spectrum was measured in a 
post VETA test where the proportional counter 
gain was set appropriately such that the high ener- 
gy continuum could be fitted. A set of apertures 
with pinhole diameters ranging from 5 |im to 20 
mm were placed at the focal plane in front of the 
XDA proportional counter in order to measure the 
encircled energy at different radii.The calculated 
encircled energy using both quick look and fitting 
techniques are presented in figure 8 . The difference 
between encircled energy calculated by each meth- 
od is very noticeable for Zr and Mo.This can be ex- 
plained by referring to the model effective area of 
the mirror in figure 6 . Notice that near Zr the effec- 
tive area curve is very steep. This implies that the 
attenuation of the continuum component in the 
XDA spectrum will be strongly dependent on ener- 
gy, resulting in a line to continuum ratio in the re- 
gion of interest quite different from the one 
measured in the BND spectrum. 

Since the largest error component in determin- 
ing the line counts in a spectrum originates from 
the continuum subtraction one desires to minimize 
the continuum counts in the ROI by selecting a ROI 
as small as possible around the X-ray line peak. Ob- 
viously for detectors with better energy resolution 
one may select narrower ROI's thus reducing the 
error due to the continuum subtraction. We have 
investigated the sensitivity of the calculated line 
events in the spectrum with the selected region of 
interest and found that it is insignificant. The varia- 
tion of the calculated encircled energy with selected 
region of interest for both the quick look analysis 
technique and the spectral fitting technique Is 
shown in figure 9. Notice that the quicklook analy- 
sis technique is very sensitive to the selected region 
of interest Initial values for the effective area used 
in the fit to the focal plane spectra were obtained by 
ray tracing the VETA-I mirror and using optical 
constants as measured by Henke et al. The next it- 
erative step is to repeat the analysis using the effec- 
tive area as determined from the fits to the spectra 
as input to the code. 


In table 2 we present the measured 1 sigma uncertainties in the XDA and BND line events. The error in 
the ratio XDA/BND at a certain energy reflects the measured uncertainty in determining the encircled en- 
ergy (EE) at that energy, since the EE for our experimental setup is approximately n x XDA/BND. 

The energy scale of the pulse height channels is also determined with an iterative procedure. We initially 
locate the pulse height channel that corresponds to the peak of the smoothed spectrum. This value is used 
as an initial guess for the fitting routine which after a certain number of iterations provides the pulse height 
channel that corresponds to the maximum of the fitted line component of the spectrum. The pulse height 












Figure 8. Measured encircled energy using both quick- 
look and fitting techniques 


treated in a separate paper by Zhao et al 6 . 


Table 2: Measured 1 sigma uncertainties in the XDA 
and BND line events 



c 

Cu 

Al 

Zr 

Mo 

°XDA 

0.37% 

0.2% 

0.17% 

0.53% 

0.68% 

°BND 

0.20% 

1.1% 

0.15% 

0.66% 

1.2% 

°X£A 

BND 

0.42% 

1.2% 

0.23% 

0.85% 

1.4% 
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Figure 9. Variation of the measured encircled energy 
\ with selected region of interest for both the quick look 
analysis technique and the spectral fitting technique 
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energy scale is then recalculated using the corrected 
value for the line peak location. This procedure is 
especially effective when dealing with low count 
spectra and spectra with large continuum compo- 
nents. 

The effects of X-ray transmission attenuation 
due to the proportional counter window mesh are 
























